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ADAPTIVE FILTER TO REDUCE MULTIPATH 

{ Cross Reference to Related Applications 
This method of reducing multipath avoids the cross term noise produced by other methods [Refs. 
1 through 10]. These other methods correlate (cross correlate or auto correlate) two versions of the 
signal each containing the direct path and multipath reflections. A correlation peak between the 
direct path and a reflection gives a measurement of the relative delay between the direct path and 
the reflected path and the amplitude and phase of the reflection. As the relative delay between the 
two signals is moved from this peak, other correlation peaks are produced when there is correlation 
between the reflections, i.e., when the relative delay equals the difference of the delays for two 
objects. This produces noise which interferes with the measurement since there are many possible 
combinations. When one delays one of the signals and subtracts it from the other with the correct 
amplitude and phase, he tends to cancel the multipath. The other reflections, however, are dragged 
along with it and produce noise. These noises are avoided by the invention described in this 

disclosure. } 

ABSTRACT 

The direct path of a radio signal from the transmitter to the receiver is frequently interfered with 
by reflections of the transmitted signal from stationary and moving objects. { This is called} These 
reflections are known as multipath noise. { This invention utilizes a new adaptive filter technique to 
reduce multipath noise. A nonlinear least squares method measures the delay, Doppler shift and 
amplitude of the multipath due to each object and subtracts a very accurate reconstruction of each 
multipath signal from the noisy signal } The multipath canceler of this invention takes as input the 
direct path signal plus multipath noise. The canceled after canceling the multipath noise, outputs a 
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purified version of the direct path signal nearly devoid of multipath noise. The output of the 
canceler is fed back to a bank of delayers and frequency shifters. The outputs of this bank of 
delayers and frequency shifters are multiplied by a set of adjustable coefficients. These coefficients 
are adjusted to form very accurate replicas of the multipath reflections caused by each reflector. 
Which replicas are subtracted from the input which is the original direct path signal plus the 
multipath noise. The result of this subtraction is called the residual and is a purified version of the 
direct path signal nearly devoid of multipath noise. The coefficients are determined by nonlinear 
least squares which finds that set of coefficients which minimize the mean square of the residual. 
Each coefficient is associated with a delay and a frequency shift. The amplitude of a coefficient is a 
measure of the magnitude of a reflector and the differential delay and frequency shift caused by said 
reflector. If an object is a target, its differential range, range rate and magnitude is got from the 
measured {multipath} differential delay, {Doppler} frequency shift and amplitude of the 
coefficient. Position and velocity of the target can be obtained by geometric triangulation with 
multiple transmitters. Target angle can be measured (by) from the relative phase angles of the 
corresponding filter coefficients , if complex, across multiple antennas receiving the same signal 
from the same transmitter. The system, comprising the antennas, receivers, signal processors, 
multipath cancelers and target estimators, can be used on a surveillance aircraft to cancel ground 
reflections and {measure} track targets. 
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ADAPTIVE FILTER TO REDUCE MULTTPATH 
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ADAPTIVE FILTER TO REDUCE MULTIPATH 
BACKGROUND 

The direct path of a radio signal from the transmitter to the receiver is frequently interfered with 
by reflections of the transmitted signal from stationary and moving objects which reflections are 
also received by the receiver. These reflected signals are delayed and frequency shifted versions of 
the direct path signal. When mixed with the direct path signal, the reflected signals corrupt the 
direct path signal and constitute what is known as {This is called) multipath noise. { This invention 
utilizes a new adaptive filter technique to reduce multipath noise. This method of reducing 
multipath avoids the cross term noise produced by other methods [Refs. 1 through 10]. } In order to 
cancel the multipath reflections, current multipath canceling methods [References 1 through 101 
first measure them. These {other} methods correlate (cross correlate or autocorrelate) two versions 
of the signal each containing the direct path and multipath reflections. For example, the direct path 
signal plus multipath reflections is multiplied by a delayed version of the signal plus multipath. The 
delay is variable. A correlation peak between the direct path and a reflection at a specific value of 
the variable delay gives a (measurement } measure of the relative delay between the direct path and 
the reflected path and a measure of the relative amplitude {and phase}of the refection. { As the 
relative delay between the two signals is moved from this peak, other correlation peaks are 
produced when there is correlation between the reflections. This produces noise which interferes 
with the measurement. When one delays one of the signals and subtracts it from the other with the 
correct amplitude and phase, he tends to cancel the multipath. The other reflections, however, are 
dragged along with it and produce noise. These noises are avoided by the invention described in this 
disclosure.} However, when the variable delay matches the relative delay between two reflections. 




undesired correlation peaks are produced. These peaks are multipath cross correlation noise and 
they corrupt the measurement process. Once the multipath reflections have been measured, the 
signal plus multipath is delayed and adjusted in amplitude such that the modified direct path signal 
portion approximately matches each reflected signal. These modified versions are, then, subtracted 
from the original signal plus multipath, reducing the multipath noise. However, even if the modified 
direct path signal portion perfectly matches the multipath reflections and they are totally canceled 
in the original signal plus multipath, the modified versions contain multipath reflections as well. 
This adds secondary multipath noise to the signal being restored and corrupts the result. The 
invention described in this disclosure avoids this multipath cross correlation noise and the 
secondary multipath noise. This invention reduces the multipath noise by an improved method of 
canceling the reflected versions of the direct path signal. This method delays a purified version of 
the direct path signal, nearly devoid of multipath noise, and cancels each component of the 
multipath noise, in contrast to other methods which delay a version of the direct path signal 
containing significant multipath noise and attempt to cancel the multipath noise with it. The 
differential delay, frequency shift, relative amplitude and phase shift due to each reflector are 
produced by the canceling process. These parameters can be used in a surveillance system to locate 
and track particular reflectors of interest, that is, targets. 
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BRIEF SUMMARY 



{ A nonlinear least squares method measures the delay, Doppler shift and amplitude of the 
multipath due to each object and subtracts a very accurate reconstruction of each multipath signal 
from the noisy signal. } This method delays a purified version of the direct path signal, nearly 
devoid of multipath noise, and cancels each component of the multipath noise, in contrast to other 
methods which delay a version of the direct path signal containing significant multipath noise and 
attempt to cancel the multipath noise with it. The multipath canceler of this invention takes as input 
the direct path signal plus multipath noise. The canceled after canceling the multipath noise, 
outputs a purified version of the direct path signal nearly devoid of multipath noise. The output of 
the canceler is fed back to a bank of delayers and frequency shifters. The outputs of this bank of 
delayers and frequency shifters are multiplied by a set of adjustable coefficients. These coefficients 
are adjusted to form very accurate replicas of the multipath noise caused by each reflector. Which 
replicas are subtracted from the input which is the original direct path signal plus the multipath 
noise. The result of this subtraction is called the residual and is a purified version of the direct path 
signal nearly devoid of multipath noise. The coefficients are determined by a method which finds 
that set of coefficients which minimize the mean square of the residual. When the residual that is, 
the signal plus multipath noise minus the replicas of the reflected signals, is minimized, the 
multipath noise is almost totally canceled. The method of determining the set of coefficients which 
minimize the mean square residual is nonlinear least squares. Newton's method for solving 
nonlinear equations was generalized by Gauss to find the parameters of planetary orbits by 
minimizing the mean of the squares of the residuals of the observations, a nonlinear problem. The 
said coefficients can be determined by the Marquardt method, an improved version of the Gauss- 
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Newton method. At first, one may think that this invention is like pulling oneself up by his own 
bootstraps. That this is not the case is proven by the mathematical explanation, the computer 
simulation in the Appendix and the fact that G. E. P. Box had very good success in a similar 
a pplication. Each coefficient is associated with a delay and a frequency shift. The amplitude of a 
coefficient is a measure of the magnitude of a reflector and the differential delay and frequency 
shift caused by said reflector. If an object is a target, its differential range, range rate and magnitude 
is got from the measured (multipath! differential delay. (Doppler) frequency shift and amplitude of 
the coefficient. Position and velocity of the target can be obtained by geometric triangulation with 
multiple transmitters. Target angle can be measured from the relative phase angles of the 
corresponding filter coefficients , if complex, across multiple antennas receiving the same signal 
from the same transmitter. The system, comprising the antennas, receivers, signal processors, 
multipath cancelers and target estimators, can be used on a surveillance aircraft to cancel ground 
reflections and {measure} track targets. 
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BRIEF DESCRIPTION OF DRAWING 



The Figure is a block diagram illustrating { the adaptive filter of} the invention to reduce 
multipath in a radio transmitted signal lj_2 { wherein the multipath is caused by reflections from 
fixed and moving objects, interferers, or targets }. The direct path of a radio signal from the 
transmitter 1 and transmitting antenna 2 to the receiving antenna 5 and receiver 6 is frequently 
interfered with by reflections of the transmitted signal from stationary 3 and moving objects 4 
which reflections are also received by the receiver 6. These reflected signals are delayed and 
frequency shifted versions of the direct path signal. When mixed with the direct path signal the 
reflected signals corrupt the direct path signal and constitute what is known as multipath noise. The 
direct path signal plus multipath noise 7, after conversion to complex samples 8, 9, 10, is whitened 
11, 12 , 13. The canceling method works best on white signals; therefore, prewhitening 12 and post 
unwhitening 18 filters are employed. The multipath canceler of this invention takes as input the 
whitened direct path signal plus multipath noise 13. The canceled after canceling the multipath 
noise 13, 14, 17, outputs a purified version 17 of the direct path signal nearly devoid of multipath 
noise. This is the residual 17. The output 17 of the canceler is fed back to a bank of delayers 11 and 
frequency shifters 16. The outputs of this bank of delayers and frequency shifters are multiplied by 
a set of adjustable coefficients 15. These coefficients are adjusted to form very accurate replicas 14 
of the multipath noise caused by each reflector. Which replicas are summed and subtracted 14 from 
the input 13 which is the whitened direct path signal plus the multipath noise. The result of this 
subtraction is called the residual 17 and is a purified version of the whitened direct path signal 
nearly devoid of multipath noise. The differential delay and {Doppler} frequency shift for each 
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{target} reflector 3, 4 is measured by the filter weight ^mn 15 which corresponds to the 

differential delay nD and (Doppler) fre quency shift caused by that {target} reflector 3, 4. 
The magnitude of the filter weight is a measure of the strength of the {target} reflector 3,4. The 

filter weight mn corresponding to a moving target 4 is a measure of the differential range 

YlDc and range rate ^fift* where C is the speed of light. The coefficients are determined by 
a method which finds that set of coefficients 12, 15 which minimize the mean square of the residual 
17. The quantity labeled residual 17 in the Figure is the received signal 13 minus the result 14 of 
application of the filter weights 12, 15 

W 2 5^3?* ' * ? %\2 ? ^13?" * * ?^21 9*^22 '^23 * # >^31 >^32 5^33?* A 
nonlinear least squares method is used to pick the filter weights which minimize the mean square of 
the residual. When the residual 17, that is, the whitened direct path signal plus multipath noise 13 
minus the replicas 14 of the reflected signals, is minimized, the multipath noise is almost totally 
canceled. The key to understanding the operation of this invention is the recognition that the 
residual 17 is a highly accurate estimate of the whitened direct path signal. This residual 17 is, then, 
unwhitened 17, 11, 18. The result is a very accurate estimate of the original direct path signal 
{ The result is a clean version of the original signal which is the received signal minus the 
multipath reflections. } { The system, with a common antenna, can be duplicated to receive two 
signals each from each of two widely separated transmitters. Geometric triangulation can be used 
to measure the two dimensional position and velocity of the target. Triplication can be used for 
three dimensional position and velocity. 



Another use of the system is to detect the angle of moving targets. If multiple antennas are 
provided, each connected with a system like that shown in the Figure, the angle of a moving target 
causing a multipath reflection to the receiving antennas can be measured by using the relative phase 
of the corresponding delay and Doppler complex coefficients across the several antennas . 

The receiving antennas, receivers and processing system can be placed in a surveillance aircraft 
the position and velocity of which is obtained by an accurate navigation system such as GPS. 
Objects on the ground, interferers and targets causing multipath reflections can be processed by the 
system. Target position and velocity can be obtained by adding the position and velocity of the 
surveillance aircraft to the measured position and velocity of the target.} 
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DETAILED DESCRIPTION 

The direct path of a radio signal from the transmitter 1 and the transmitting antenna 2 to the 
receiving antenna 5 and receiver 6 is frequently interfered with by reflections of the transmitted 
signal from stationary 3 and moving objects 4 which reflections are also received by the receiver 6. 
These reflected signals are delayed and frequency, or Doppler, shifted versions of the direct path 
signal. When mixed with the direct path signal, the reflected signals corrupt the direct path signal 
and constitute what is known as multipath noise. 

The direct path signal plus multipath noise 7, after conversion to complex samples 8, 9, 10, is 
whitened 11, 12 , 13. The multipath canceler of this invention takes as input the whitened direct 
path signal plus multipath noise 13. The canceler, after canceling the multipath noise 13, 14, 17, 
outputs a purified version 17 of the direct path signal nearly devoid of multipath noise. This is the 
residual 17. The output 17 of the canceler is fed back to a bank of delayers 11 and frequency 
shifters 16. The outputs of this bank of delayers and frequency shifters are multiplied by a set of 
adjustable coefficients 15. These coefficients are adjusted to form very accurate replicas 14 of the 
multipath noise caused by each reflector. Which replicas are summed and subtracted 14 from the 
input 13 which is the whitened direct path signal plus the multipath noise. The result of this 
subtraction is called the residual 17 and is a purified version of the whitened direct path signal 
nearly devoid of multipath noise. The coefficients are determined by a method which finds that set 
of coefficients 12, 15 which minimize the mean square of the residual 17. When the residual 17, 
that is, the whitened direct path signal plus multipath noise 13 minus the replicas 14 of the reflected 
signals, is minimized, the multipath noise is almost totally canceled. The key to understanding the 
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operation of this invention is the recognition that the residual 17 is a highly accurate estimate of the 
whitened direct path signal. This residual 17 is, then, unwhitened 17, 11, 18. The result is a very 
accurate estimate of the original direct path signal. 

The method of determining the set of coefficients 12, 15 which minimize the mean square 
residual 17 is nonlinear least squares. Newton's method for solving nonlinear equations was 
generalized by Gauss to find the parameters of planetary orbits by minimizing the mean of the 
squares of the residuals of the observations, a nonlinear problem. The said coefficients can be 
determined by the Marquardt method, an improved version of the Gauss-Newton method of 
nonlinear least squares. At first, one may think that this invention is like pulling oneself up by his 
own bootstraps. That this is not the case is proven by the mathematical explanation, the computer 
simulation in the Appendix and the fact that G. E. P. Box had very good success in a similar 
application. 

Each coefficient 15 is associated with a delay and a frequency shift. The amplitude of a 
coefficient is a measure of the magnitude of a reflector and the differential delay and frequency 
shift caused by said reflector. If an object is a target, its differential range, range rate and magnitude 
is got from the measured differential delay, frequency shift and amplitude of the coefficient. 
Position and velocity of the target can be obtained by geometric triangulation with multiple 
transmitters. Target angle can be measured from the relative phase angles of the corresponding 
complex filter coefficients across multiple antennas receiving the same signal from the same 
transmitter. The system, comprising the antennas, receivers, signal processors, multipath cancelers 
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and target estimators, can be used on a surveillance aircraft to cancel ground reflections and track 
targets. 



In order to cancel the multipath reflections, current multipath canceling methods [References 1 
through 101 first measure them. To do the measurement, { This method of reducing multipath 
avoids the cross term noise produced by other methods [Refs. 1 through 10]. }these {other} 
methods correlate (cross correlate or autocorrelate) two versions of the signal each containing the 
direct path and multipath reflections. For example, the direct path signal plus multipath reflections 
is multiplied by a delayed version of the signal plus multipath. The delay is variable. A correlation 
peak between the direct path and a reflection at a specific value of the variable delay gives a 
{measurement} measure of the relative delay between the direct path and the reflected path and a 
measure of the relative amplitude {and phase} of the refection. However, when the variable delay 
matches the relative delay between two reflections, undesired correlation peaks are produced. 
These peaks are multipath cross correlation noise and they corrupt the measurement process. 
{ As the relative delay between the two signals is moved from this peak, other correlation peaks are 
produced when there is correlation between the reflections. This produces noise which interferes 
with the measurement. When one delays one of the signals and subtracts it from the other with the 
correct amplitude and phase, he tends to cancel the multipath. The other reflections, however, are 
dragged along with it and produce noise. These correlation noises are avoided by the invention 
described in this disclosure.} 

Once the multipath reflections have been measured, the signal plus multipath is delayed and 
adjusted in amplitude such that the modified direct path signal portion approximately matches each 
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reflected signal. These modified versions are, then, subtracted from the original signal plus 
multipath, reducing the multipath noise. However, even if the modified direct path signal portion 
perfectly matches the multipath reflections and they are totally canceled in the original signal plus 
multipath, the modified versions contain multipath reflections as well. This adds secondary 
multipath noise to the signal being restored and corrupts the result. 

The invention described in this disclosure avoids this multipath cross correlation noise and the 
secondary multipath noise. This invention reduces the multipath noise by an improved method of 
canceling the reflected versions of the direct path signal. This method delays a purified version of 
the direct path signal, nearly devoid of multipath noise, and cancels each component of the 
multipath noise, in contrast to other methods which delay a version of the direct path signal 
containing significant multipath noise and attempt to cancel the multipath noise with it. The 
differential delay, frequency shift, relative amplitude and phase shift due to each reflector are 
produced by the canceling process. These parameters can be used in a surveillance system to locate 
and track particular reflectors of interest, that is, targets. 

The Figure is a block diagram illustrating the adaptive filter of the invention to reduce multipath 
in a radio transmitted signal 1^2 wherein the multipath is caused by reflections from fixed 3 and 
moving objects 4, interferers, or targets 4. 

{ The method works best on white signals; therefore, prewhitening and post unwhitening filters are 
employed. 
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The differential delay and Doppler shift for each target is measured by the filter weight which 
corresponds to the differential delay nD and Doppler shift mf caused by that target. The magnitude 
of the filter weight is a measure of the strength of the target. 

The quantity labeled "Residual" in the Figure is the received signal minus the result of application 
of the filter weights w x , w 2 , w 3 , jc 1 , , x n , jc 13 • • • , x 2} ,x 22 ,x 2 y x 3} , x 32 , x 33 • • • .A nonlinear least squares 
method is used to pick the filter weights which minimize the mean squares of the "Residual". The 
result is a clean version of the original signal which is the received signal minus the multipath 
reflections.} 

{Refer to the Figure} The original signal is modulated and amplified by the transmitter, Tx 1 
[Ref 17, ROHDE, p. 14]. The original signal is defined as that signal which rides on the carrier. If 
one band passes the transmitted signal and, then, band shifts the result down near zero frequency, 
he obtains the original signal. The original signal is formed after the unique modulation such as 
AM, FM, etc. is done. After transmission from the antenna 2, most of the power takes the direct 
path to the receive antenna 5. Some power, however, is reflected from stationary 3 and moving 4 
objects such as buildings and aircraft. When a reflection arrives at the receive antenna 5, it is a 
delayed and attenuated version of the direct path signal. If the reflector is moving 4, there is also a 
frequency shift. At the receive antenna 5, these reflections combine with the direct path and cause 
{This tends to add or subtract from the direct path causing } what is known as multipath 
interference noisefRef. 18, SKOLNflC,p. 18.46]. 
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This invention removes the multipath noise from the received signal. The direct path signal plus 
the multipath noise plus the receiver noise, after conversion to complex samples 10, is whitened in 
eq. (15) 13. The multipath 14 is canceled 13, 14, 17 in eq. (17) by nonlinear least squares in eq. 
(20). Unwhiteninfi 17, 11, 18 is performed in eq. (16) to produce the restored direct path signal. The 
processing which accomplishes this will be described in detail and the operation will be explained 
in what follows. 

The signal at the receive antenna 5 is amplified and demodulated by the receiver, Rx 6 [Ref 17, 
ROHDE, p. 361 ][Ref. 19, VAN VALKENBURG, p. 23-6]. This demodulation consists of band 
passing the received signal and band shifting, hetrodyning, the result down to neiar zero frequency. 
At this point 7, the direct path signal plus multipath noise plus receiver noise is analog. It is, then, 
fed to Analog to Digital processing 8. Which processing is comprised of a presampling filter, 
sampler and analog to digital converter { It is then passed through a presampling filter, sampled and 
converted from analog to digital } [Ref. 19, VAN VALKENBURG, p. 28-7]. The output of Analog 
to Digital processing is a sequence of real digital samples which are fed to Real to Complex 
processing 9. An alternative to the analog receiver 6 is a digital receiver. In which case, the Analog 
to Digital processing is moved forward toward the antenna 5. 

The preferred initial embodiment of this invention is to process the output of the Analog to 
Digital conversion with a digital computer and to produce the final results in digital form. The 
restored signal will be in digital form which can be converted to analog if desired. Later 
embodiments would include special purpose signal processing hardware for the signal processing in 
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the Figure: Real to Complex, Whitening, Multipath Canceler and Restore Spectrum. The nonlinear 
least square function could also be done in special purpose hardware. 

Real to Complex processing 9 converts the sequence of real samples to a sequence of complex 
samples. The {These} real digital samples have a Fourier transform the real part of which is 
symmetric about zero frequency and the imaginary part of which is antisymmetric about zero 
frequency. From this redundancy, one can see that one needs only half of the frequencies in the 
Fourier transform, that is, half the time samples if they are complex. The real samples have a 

maximum frequency of IF and a sampling frequency of 4F . The complex samples will have 

__{_ 

a maximum frequency of r and a sampling frequency of where D is the 

time between complex samples. To convert from real samples to complex samples 9, the real 

D 

-jlTlFt t - l 

samples are first frequency down shifted by multiplying by ^ where { I } 9 



{ F = yd and y is the sampling frequency } is time[Ref 14, PAPOULIS, p. 15] and / is the real 
sample index. These complex samples are then passed through a low pass digital filter with a 
maximum frequency of F [Ref. 16, RABINER, p. 100][Ref. 13, OPPENHEIM, p. 261]. The 

1 

2F = — 

output of the low pass filter is resampled at a frequency of { y D , i.e., } , that is, 

every other sample is taken to obtain the complex samples of the received signal 10. 
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The reason for converting to complex samples 10 is that, during processing, the samples will not 
only be delayed 11 but also frequency shifted 16. This frequency shifting can be easily affected by 
multiplying the complex samples by a complex exponential 16. See eqs. (L 12, 13, 17). Complex 
phase is also used in measuring the angle of a target See eq. (38). 
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The received signal, that is, the direct path signal plus multipath noise plus white receiver noise, 
after demodulation, analog to digital conversion and conversion to complex samples 10, can be 
modeled as 

N M 
n-\ m-- M 

where 



$i = original direct path signal (2) 

I = sample index (3) 

D = sampling time interval (4) 

Yl = time delay index for each reflector (5) 

f = {Doppler} frequency shift frequency interval (6) 

YYl = {Doppler} frequency shift index for each reflector (7) 

^mn = Multipath amplitude coefficient for each reflector (8) 

J = unit imaginary number (9) 

@ = base of natural logarithms (10) 

fl"- P i £11} 

E = e il4 ' 021 
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t = ID £13} 

^li = white receiver noise (14) 
Each term in the double sum of eg. (1) represents a reflector causing a delay of the direct path 

signal of , a frequency shift of 

mf 

and an amplitude attenuation of mn . When 

E of eq, (12) and t of eq. (13) are substituted into eq. (1), one obtains the term in the double 

n c p j27nnf{i-n)D _ , 

sum_ /w« /-« which represents the I *™ multipath sample at the time 

delay of the direct path signal of nD and the frequency shift of 

mf 

. The complex 

exponential coefficient accomplishes the frequency shift of the time delayed complex direct path 

signal ^i-n • The coefficient ^mn is the complex magnitude (at the receive antenna) of the 
multipath reflector at differential (between the direct path length and the reflected path length) 

distance ftPc and differential velocity Wife w here C is the speed of light 

If the spectrum of the direct path signal is not white, there can be correlation between one time 
point of the signal and another. This allows the multipath canceled with the appropriate delay 
coefficients 15, to partially cancel the direct path signal as well as the multipath. Putting the direct 
path signal through a whitening filter 12 avoids this signal suppression. This will be explained in 
eq. (26). Note that since the multipath reflections are replicas of the direct path signal, they are also 
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whitened.. Note that the initial white receiver noise is modified by the whitening filter and is no 



longer white. 



The direct path signal plus multipath plus white receiver noise 



(7/ 10 is fed into the whitening 



filter. The whitening filter is comprised of a sequence of one sample delays 11. This sequence is 
ta pped after each delay and fed to the filter weights ^/ 12. These inputs are multiplied by the 



These filter weights are adjusted by least squares to remove any correlation from the input 10. 
This adjustment is done as part of the nonlinear least square minimization of the residuals. The 
identical (to the whitening filter) filter weights 18 are used in the unwhitening filter to restore the 
spectrum. The operation of the whitening filter will be explained in eq. (23). 

The output of the whitening filter 13 is fed to the multipath canceler which subtracts off the sum 
of the multipath replicas 14 to form the residual 17. The key to understanding the operation of this 
invention is the recognition that the residual 17 is a highly accurate estimate of the whitened direct 
path signal. The residual is passed through the unwhitening filter 11, 18 which restores the 



filter weights the outputs from which are summed and subtracted from (7/ 10 forming the 



whitened direct path signal plus multipath plus receiver noise 13 



L 
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spectrum of the direct path signal yielding a highly accurate estimate of the original direct path 



The unwhitening filter is similar to the whitening filter, however, it is a recursive filter, that is, the 
outputs are fed back as inputs. The sequence of delays 11 is the same and the filter weights 18 are 
identical to those of the whitening filter 12. The operation of the unwhitening filter will be 
explained in eqs. (31) and (32). 

The replicas 14 are formed as follows. The output of the multipath canceler, the residual 17, is a 
very good estimate of the whitened direct path signal. It is fed back to a bank of delayers 11 and 
frequency shifters 16. The outputs of the bank of delayers and frequency shifters are multiplied by a 
set of adjustable coefficients 15. These coefficients are adjusted, along with the whitening 
coefficients, to form the very accurate replicas of the whitened multipath noise caused by each 
reflector. These are summed and subtracted 14 from the multipath canceler input 13 to form the 
residual 17. The coefficients 

W| ? 9 ? • • • 9 j 9 ? 9 # • # ? 9 %22 * *^23 > " * # * *^31 * ^32 9 ^33 9 * * are 

determined by minimizing the root mean square (r.m.s.) residual 17 by the Marquardt method of 
nonlinear least squares. The minimization of the residual 17 by nonlinear least squares will be 
described in eq. (20). The operation of the multipath canceler will be explained in eqs. (24-29). 



signal 



L 
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The differential delay and frequency shift for each target is measured by the filter weight ^mn 

which corresponds to the differential delay nD and frequency shift caused by that target. 
The magnitude of the filter weight is a measure of the strength of the target. 

From the multipath canceler processing shown in the Figure, the residual 17 is 

N M 

n = V /-Z ;/ = l,2,-,/ Q7} 

n=\ m=-M 



where v i 13 is the whitened direct path signal plus multipath plus receiver noise and the mn 



15 are the frequency shift-delay filter coefficients. The last term in eq. (17) is the sum of the 
multipath replicas 14. Note that the residual is defined recursively. To start the recursion, one 

*• • 1 L N M 

assumes that r i- n is zero for I ~ W < 1 L { r t .= q i , - X - X X x mn r i-n Em 

1=1 n=\ m=~M 

where the w t are whitening filter coefficients and the are the delay-Doppler filter coefficients. 
This equation } Equations (15) and (17) correspond to that of Box [Ref 1 1, BOX, p. 498]. One can 
minimize the root mean squared (r.m.s.) residual 



1M 
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over 

W l > W 2 > ' # # > > X - A/,1 ? X -M,2 5 * * * ? X -M,N 9 X - > 

Y • • • Y • • • Y 

A -M+l,2> > A -A/+1,A" ^^M 9 N 

by the Marquardt method of nonlinear least squares [Ref. 11, BOX, p. 504, p. 505] [Ref. 12, 
MARQUARDT, p. 431], 
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The Marquardt method is a sophisticated form of simple nonlinear least squares. Simple 



nonlinear least squares can be expressed 



x f = x- 



a r 1 
\d xj 



or' 
\d x) 



-i 



a r 



\d x) 



I2Q} 



where the 



l(L + N{2M + \)) A 



dimensional coefficient vector is 



x = Of ,wf ,w* ,w J 2 , — 5 <^-Vl^-M,P 



x R x J -~x R x J x R x J 

*-M,2 ' A -M,2 9 -M,N ? -M,N ' ' — Af+1,1 > 



*-M+1.2 9 X -M+1.2 ?" ' * 9*-M+lJV 



9-* / -M+\,N9"°- x M,N9- x M,N ) (21) 



2/ 



the £1 dimensional residual vector is 



r = (r 1 *,r 1 J ,r 2 R ,r 2 J ,-,rj l ,r/)\ 



(22) 



and the values of ^/ are found from eqs. (15) and (17) using X 12, 15 and 10, the sampled 



complex direct path signal plus multipath noise plus white receiver noise. The superscript 



R 



means "real part" and the superscript J means "imaginary part." The dimensionality of Y must 



be greater than that 



i fJC/rhe2/b^2(Z + M2M + l)) 



partial derivative matrix 
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is calculated numerically, the superscript 



means matrix transpose, the superscript 



-1 



means matrix inverse and X is the new value of X in the iterative nonlinear least 



squares calculation of eg. (20). 

The method of nonlinear least squares is used to find the values of those coefficients which 
minimize the mean square residual. The method is iterative, that is, one starts with a set of 
coefficients, calculates the residuals, calculates an increment for each coefficient and increments 
the previous set of coefficients to find a new set of coefficients. See eg. (20). This process is 
repeated until it converges, that is, until one gets the same answer twice. 

For example, let the coefficients be 

^1 9 ^2 ?^12 >^13 5^21 5*^22 5*^23 '^31 ?^32 ?^33 , where there are 3 

delays in the whitening filter and 3 delays and 3 frequency shifts in the canceler. There are a total of 
12 coefficients to be determined, and since each complex coefficient has a real and an imaginary 
part, a total of 24 values to be determined. The input signal plus multipath noise has been converted 
to a series of complex digital samples, say 1000 of them, that is, 2000 values. One starts the process 
by choosing a set of coefficients, say, all zero. The sequence of input samples are whitened under 

the control of the W coefficients and the multipath replicas 14 are subtracted from them to form 
the sequence of 1000 residual samples. The sequence of residual samples is fed back into the bank 
of delayers and frequency shifters controlled by the X _ coefficients. Since the residual is an 
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approximation to the whitened direct path signal the outputs of the bank of delayers and frequency 
shifters are sequences of whitened multipath replicas. These are subtracted from the whitened direct 
path signal plus whitened multipath. The result of this process is a sequence of 1000 complex 
residuals, that is, 2000 values. For the first iteration, the multipath may not be reduced, if at all. 
However, for the second iteration, an increment of the coefficients is calculated by nonlinear least 
squares. The first step is to calculate the partial derivatives of the 1000 residuals with respect to the 
12 coefficients. This is done numerically by stepping each coefficient, one at a time, by a small 
delta and recalculating the 1000 residuals for each coefficient. The partial derivative is found by 
taking the difference between the stepped residual and the unstepped residual and dividing this 
difference by the delta. This gives a partial derivative matrix of 1000 by 12 complex elements, that 
is, 2000 by 24 real elements. A square 24 by 24 matrix is formed by matrix multiplying the 24 by 
2000 transpose of the partial derivative matrix by the 2000 by 24 partial derivative matrix itself 
The 24 increments to the first set of coefficients are, then, found by matrix multiplying the 24 by 24 
matrix inverse of the square matrix times the 24 by 2000 transpose of the partial derivative matrix 
times the 2000 element vector of residuals. This process is repeated, that is, iterated until it 
converges, that is, until the same result is got twice. Upon convergence, the multipath content of the 
input signal plus multipath has been canceled and the residual, after the spectrum has been restored, 
is a purified version of the direct path signal. 
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The following is a mathematical explanation of how the invention removes the multipath and 
restores the original signal. First consider the whitening process. For demonstration purposes, let 

all the ^mn coefficients 15 in eq. (17) be zero. Then, the output of the whitener 13 equals 
the residual 'j" 17. When one solves eq. (15) for (7/ 10, the direct path signal plus multipath 



noise plus white receiver noise, he gets 
L 

am n 4- r 

(22) 



1=1 



The process minimizes the root mean square of 17 over ^/ 12 . The result is an optimum 

linear predictor of Qj , as the current value is predicted by a linear combination of previous 
Y 

values. The residual ' / is the prediction error and must be white, for if it were not, it would 
contain some correlation which could be used to further reduce the prediction error. Now, since 

V Y 

i is the output of the whitener and is equal to ' / , the output 13 is white. 

{ To show that this removes the Multipath and restores the original signal } 

Now, to explain multipath canceling, using eq. (15), define the whitened direct path signal as 



D . = S — / W,S ; 

Fi z I i-l (24) 

and the modified white receiver noise as 
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Then, substituting eq. (1) into (15), using eqs. (24) and (25), and substituting eq. (15) into eq. (17), 
one obtains the form of the residual {is} 

N M N M 

n=l m=-M n~l m=-M 

126) 

{ The root mean square residual is minimized when 

x rm ^o mn e- J2mnfhD 
and 

Note that only delayed samples of the whitened direct path signal Pi {, i.e., p i _ n 5 n * 0, }, 
contained by ^i-n * are subtracted from Pi , so that it is not possible to cancel the whitened 

direct path signal Pi by varying ^mn since Pi , because it is white , is uncorrected with 

delayed versions of itself. Note that ^*i-n contains Pi-n but not A, 

{ Adding the above two equations which describe the action of the whitening filter, one obtains 

L 

n «(*i + %)-Z w /( J /-/ + '7w) 

Solving this equation, one gets 
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#1 



/=1 



which, from the Figure, gives the } 
In eq. (26) if 

X mn ~ a mn e (27) 



and 

r i = Pi (28) 

that is, if the modified receiver noise ^Di is zero, the fourth term cancels the second term exactly 
and the multipath is totally removed. This occurs when the root mean square of the residual _is 
minimized. This is so because Pi is white and cannot be reduced further by linear combinations 
of ^*i—n since ^i-n contains Pi-n but not Pi . Of course, the modified receiver noise 

is small, but not zero. Eqs. (27) and (28) are only approximate and the multipath is not totally 
canceled. This cancellation error ^ , the difference of the two double sums in eq. (26), due to a 
small , is small and eq. (26) becomes 

r, =P, + 6+5 (291 
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Note that eg. (29) shows that the residual is a good approximation to the whitened direct path 

signal Pi . The key to understanding the operation of this invention is the recognition that the 
residual is a highly accurate estimate of the whitened direct path signal. 

From the definition of unwhitening in eq. (16), define the unwhitened cancellation error 



$ = $ + Z w i s *-i (30) 



Solving (24) and (25) for and ^7/ , and adding ^ ^7/ and &j . one gets 

L 

s i + fit +£,= P, + ft + 5 + Z w i( s t-i + Vi-i + £ i-i) (3n 

Comparing eq. (3 1) to eq. (16) and using eq. (29), one has, for the restored direct path signal, 

u i = s, + rii + e, (32) 

which is the output of the unwhitening filter. This is the desired result which restores the original 

direct path signal without multipath and with a small error ^7/ . The quantity js 

the original signal in complex form. To obtain the unmodulated original signal, one must, then, 
demodulate the unique modulation, such as AM, FM, etc., which occurred before the original signal 
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was formed. If a reflector is a target the ^mn imply target differential range, range rate and 
magnitude. 
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The Appendix contains listings of computer programs written in the APL language and a sample 
computer run. The programs embody the basic features of this invention. The computer run is a 
demonstration of its performance on simulated real data. 

The program NORM generates simulated data. The program DELRES calculates the delay 
residuals. The program NLS minimizes the residuals defined in DELRES by nonlinear least 
squares. 

The program NLS is first edited so as to output the residuals as RR. The program NORM is used 
to generate 100000 mean zero, variance L real, normally distributed, random numbers in X5. This 
represents a white 50000 Hz broadband signal sampled at 100000 samples per second with r.m.s. 
(root mean square) of 1. Two multipaths are simulated both of 0. 1 magnitude. One is delayed 1 
sample; the other is delayed 3 samples and both are added to X5. The result is put into XX5. The 
signal to multipath noise ratio is calculated to be 17 dB. The NLS program is applied to XX5 and 
the two multipath noises are canceled to a level of 51 dB signal to multipath noise ratio. 

A simple test of {this} the concept of this invention was made with the computer program 
shown in the Appendix. The Marquardt method, similar to eq. (20), for nonlinear least squares [Ref 
11, BOX, p. 504, p. 505] [Ref. 12, MARQUARDT, p. 431] was programmed in the APL language 
in program NLS. The equation for the residual, eqs. (15) and (17), was programmed in DELRES 
with 




when 



m*0 



(33} 



and 
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{If 

w, *0 

this corresponds to the equation of Box [Ref. 1 1, BOX, p. 498] and is a general recursive digital 
filter [Ref. 15, PAPOULIS, p. 45].} 

That is, the reflectors were assumed to be stationary, causing no frequency shift, and the direct path 
signal was assumed to be white and did not need whitening. 



The original signal ^/ was generated with a Gaussian , or normally distributed, random number 

generator in program NORM with root mean square (r.m.s.) of 1 and 7 = 100000 .This 
corresponds to one second of a 50000 Hz bandwidth signal represented by 100000 real samples. 
The data is white, zero mean, Gaussian with r.m.s. = 1. 



The received signal (7/ was generate d, eq. (1), with 

^01 = ^03 = , otherwise = 0 (35) 
and 

This corresponds to two multipaths both of magnitude 0. 1 , one is delayed 1 sample, the other is 
delayed 3 samples. The white signal itself provides an effective noise which perturbs the estimates. 



34 



fl 



The process converged in three iterations, the final r.m.s. residual = 1.002370568. The three 
coefficients ^01 ' ^02 ' *^03 are shown to be approximately = 0.1, 0.0, 0.1 matching 

^01 > *^12 > ^03 with r.m.s. errors of about 0.003. The mean residual was about 0.0087. The 
original signal to multipath noise ratio was 17 dB. After processing, the two multipaths were 
canceled to a signal to multipath noise ratio of 51 dB. 

The filter weight ^mn corresponding to a moving target is a measure of the differential range 



antenna, can be duplicated to receive two signals each from each of two widely separated 
transmitters. Geometric triangulation [Ref. 18, SKOLNIK, p.25.5, p.25. 13] can be used to measure 
the two dimensional position and velocity of the target. Triplication can be used for three 
dimensional position and velocity. 

Given the differential path lengths from two transmitters to a target to the receiver, one can 
calculate the two dimensional position of the target from geometry. The positions of the 
transmitters and the receiver are known. The differential delays, therefore, are functions of the two 
dimensional target position. For example, the length of the direct path from one of the transmitters 
to the receiver can be calculated from their known positions. The length of the reflected path 
depends on the target position. One adds the distance from the transmitter to the target to the 
distance from the target to the receiver. The differential path length is the reflected path minus the 
direct path. There is also a measurement of the differential path length. One, therefore, has two 




where C is the speed of light. The system, with a common 
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equations, one for each of the two measurements, in two unknowns, the two dimensions of the 
position of the target. These two nonlinear equations can be solved for target position. The three 
dimensional position of the target can be found with three measurements from three transmitters, 
that is, by solving three equations in three unknowns. Target velocity can be obtained in a similar 
manner from the differential range rate measurements. 

Another use of the system is to detect the angle of moving targets. If multiple receive antennas 
are provided, each connected with a system like that shown in the Figure, the angle of a moving 
target causing a multipath reflection from a single transmitter to the target to the receiving antennas 

can be measured by using the relative phases of the corresponding differential delay nD and 

fflf x 

{Doppler} frequency shift "u complex coefficients ^mn _15 across the several receiving 

antennas [Ref. 18, SKOLNIK, p.3.34, p.3.35]. 

For example, the differential complex phase angle at a receive antenna is 

>=f(* + 4.-*0 cm 



where ^ is the wave length of the transmission, R is the distance from the transmitter to the 

R R 

target 1 is the distance from the transmitter to the receive antenna and 2 is the distance 
from the target to the receive antenna. The difference in the phase angle between this receive 
antenna and another antenna located at a distance of CI away is 
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08} 



where ^1 is the angle from the line between the two receive antennas to the transmitter and 



2 * the desired angle to the target. 



The receiving antennas, receivers and processing system can be placed in a surveillance aircraft 
the position and velocity of which is obtained by an accurate navigation system such as GPS. 
Objects on the ground, interferers and targets causing multipath reflections can be processed by the 
system. Ground reflections, known as clutter, can be canceled by this invention as, from the moving 
platform, they appear to be moving. The position and velocity of a moving target can also be 
measured with respect to the moving platform. Target position and velocity with respect to the 
ground can be obtained by adding the position and velocity of the surveillance aircraft to the 
measured position and velocity of the target. 



a 



2 is the angle from the line between the two receive antennas to the target. Since CI and 



a 



1 can be calculated from known geometry and 



is measured, eq. (38) can be solved for 
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REFERENCE TO APPENDIX CONTAINING COMPUTER PROGRAM 



The Appendix contains listings of computer programs written in the APL language and a sample 
computer run. The programs embody the basic features of this invention. The computer run is a 
demonstration of its performance on simulated real data. 

The program NORM generates simulated data. The program DELRES calculates the delay 
residuals. The program NLS minimizes the residuals defined in DELRES by nonlinear least 
squares. 

The program NLS is first edited so as to output the residuals as RR. The program NORM is used 
to generate 100000 mean zero, variance 1, normally distributed, real random numbers in X5. This 
represents a white 50000 Hz broadband signal sampled at 100000 samples per second with r.m.s. 
(root mean square) of 1 . Two multipaths are simulated both of 0. 1 magnitude. One is delayed 1 
sample; the other is delayed 3 samples and both are added to X5. The result is put into XX5. The 
signal to multipath noise ratio is calculated to be 17 dB. The NLS program is applied to XX5 and 
the two multipath noises are canceled to a level of 51 dB signal to multipath noise ratio. 
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APPENDIX 



COMPUTER PROGRAMS AND SAMPLE RUN 



VN0RMIU1V 

[0] X+NORM NiZiTiR 

[1] a 

[2] A 961204.1712 

[3] p GENERATE N NORMAL RANDOM VARIABLES, MEAN 0, SIGMA 1 

[4] p URL IS SEED 
[5] Z«-2,l"tf*2 

[6] Z+Zp( *10000 )x? ( x/Z)pl0000 

[7] r«-(02 )xZ[l; ] 

[8] R*-( -2x®Z[2; ] )*0.5 

[9] jsr-«-jir+ , ( i io.xi?)x2 lo.or 

V 1996-12-10 11.00.01 ( GMT-h ) 

VZ?£Li?£S[[]]V 

[0] J?«-r DELRES X\N\I 

[I] R ;; 

[2] p 000323.2241 

[3] R DELAY RESIDUALS 
[4] JIH-p* 
[5] R+0*Y 

[6] flCiJin+rcijn 

C 7 ] J-«-0 

[ 8 ] NEXTI : I+I+ 1 

[9] RLN+I]+YlN+n-RLN+I-iNl + .xX 
[10] +NEXTIx\I<(.fiY)-N 

[II] J?-*-iVr+ J? 

7 2000-03-25 22.51.47 (GMT-*) 

VJ7L5[D]V 

[0] AT^ZU RES;X',bXiR;PiJibR;AiB;D;EMiEiI;EE',AA 

111 fl ; ; 

[2] R 800225.0914 800820 000321.1532 

[3] fl NONLINEAR LEAST SQUARES 

[4] r CHAR. VECT. RES IS RESIDUAL FUNCTION NAME 

[5] X+XDlliI 

[6] LX+XD12;] 

[ 7 ] ' 1 

[8] 1+0 

[9] NBXTI:I+I+1 

[10] R+tRES, • AT' 

[11] fl (J? A B)+tRES,' (X LX ) • 

[12] fl (*/ + /(i?*2),[0.5]l)*0. 5 

[13] (*/ + /(J?*2),[0.5]l)*0.5 

[14] fl ANALYTIC RESIDUAL PARTI ALS 

[15] fl P+t{RES, »P» ), • AT' 

[16] n-*- AMP AS 
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[17] A CALCULATE R, A AND B 

[18] A (J? A B)+*RES t ' (X A*)' 

[19] A +CALCRAB 

[20] A NUMERIC RESIDUAL PARTI ALS 

[21] P«- • ' 

[22] J+0 

[23] NEXTJ-.J+J+l 

[24] bR+(*RES t > X+hXtJlxJ=xpX' )-R 

[25] P«-P,Ai?*A;r[J] 

[26] +NEXTJx\J<pX 

[27] P<-!$K(p;n,(pP)*p*)pP 

[28] 4J7P;li?:>l«-($P)+.xP 

[29] fl«-($P)t.xJ? 

[30] CALCRAB: 

[31] Z?«-*(((0 l + p4)p4)[ ; 1]+1£"10 )*0 . 5 

[32] AA+(AxDo . xZ? ) + 1^-10 x ( P 4)p 1,0x4 

[33] ;r-eAr-Px(2?xfl)iiiii 

[34] ^jrwrjxijo 

[35] R+tRES, ' 

[36] A (R A B)«-*RES, ' (X AAT)» 

[37] EM++/+/R, [0.5]1 

[38] £«-(*/t/(J?*2),[0.5]l)*0.5 

[39] £M),[0. 5 ]£x(((0 l + p>lA ) p ( D° . xD ) xftAA ) [ ; 1 ] , 1 ) *0 . 5 

[40] A EE,E 

[41] E 

[42] A ' » 

[43] A XE 

[44] A +NEXTIxiI<20 

[45] A RR+R 

V 2000-03-29 23.00.55 (GWr-4) 

VJITL5[D45] 

[45] A i?J?-«-J? 

[45] RR+R 

[46] V 



p*5«-.flr0j?M 10 00 00 

100000 

OS«-(*/ + /((AT5)*2),[.5]l)*.5 
1.002357739 

pZZ5-<-Ar5[3+i99997] + .lxAr5[( X99997 )o . + -l-n4] + .xl 0 10 

99997 

[>#«-(*/ + /( (XX5-3+X5 )*2),[.5]1)*.5 
0.1418755036 

20xlO®S*J7 
16 . 98230664 

(2 3p0 0 0 .001 . 001 .00DNLS ' XX5 DELRES' 
1.011856929 
1.002412348 
1.00237057 
1.002370568 
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0.09896411424 0.001084301862 0.09780701046 0.008745036907 
0.003147221067 0.003163105505 0.003146911359 1.002370568 

□ (*/ + /(( RR - 6 + X 5 )* 2 ), C . 5 ] 1 )*. 5 
0.002841730873 

20xlO®S*tf 
50.94879605 
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CLAIMS 



What is claimed is: 

8. A communication system comprising a radio transmitter transmitting a signal which is interfered 
with by fixed reflectors and moving reflectors, said system comprising an antenna, a receiver and a 
multipath reduction subsystem comprising an analog to digital converter, a real to complex 
converter, a whitening filter with complex coefficients and a multipath canceller comprising an 
array of complex delay-Doppler shift coefficients, delayers of one sample and shifters of one 
frequency increment, said array producing a residual from which the interferences from said fixed 
reflectors and said moving reflectors have been removed, said multipath canceller comprising a 
minimizer minimizing the mean square of said residual over the said whitening filter coefficients 
and the said multipath canceller delay-Doppler coefficients, said subsystem comprising a spectrum 
restoring filter utilizing the said whitening filter coefficients restoring said original signal. 

9. A multipath reduction subsystem for use in a communication system comprising a radio 
transmitter transmitting a signal which is interfered with by fixed reflectors and moving reflectors, 
said system comprising an antenna, a receiver and said multipath reduction subsystem comprising 
an analog to digital converter, a real to complex converter, a whitening filter with complex 
coefficients and a multipath canceller comprising an array of complex delay-Doppler shift 
coefficients, delayers of one sample and shifters of one frequency increment, said array producing a 
residual from which the interferences from said fixed reflectors and said moving reflectors have 
been removed, said multipath canceller comprising a minimizer minimizing the mean square of 
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said residual over the said whitening filter coefficients and the said multipath canceller delay- 
Doppler coefficients, said subsystem comprising a spectrum restoring filter utilizing the said 
whitening filter coefficients restoring said original signal. 

10. A multipath reduction method for use in a communication system comprising a radio 
transmitter transmitting a signal which is interfered with by fixed reflectors and moving reflectors, 
said system comprising an antenna, a receiver and an analog to digital converter, said multipath 
reduction method comprising reducing multipath by converting real sampled data to complex 
sampled data, whitening said complex sampled data by filtering with complex coefficients and 
cancelling said multipath by delaying and frequency shifting the residual under control of the 
complex delay-Doppler shift coefficients producing said residual from which the interferences from 
said fixed reflectors and said moving reflectors have been removed, said whitening filtering 
coefficients and said multipath cancelling delay-Doppler shift coefficients being determined by 
minimizing the mean square of said residual, said multipath reduction method comprising restoring 
the spectrum of said signal by filtering utilizing the said whitening filtering coefficients to restore 
said signal. 

1 1 . A surveillance system with targets which are moving reflectors, with interfering fixed reflectors 
and with interfering moving reflectors, said system comprising one or more stationary or moving 
transmitters, one or more stationary or moving antennas, one or more receivers and a surveillance 
subsystem comprising a multipath reduction subsystem for each receiver, said multipath reduction 
subsystem comprising an analog to digital converter, a real to complex converter, a whitening filter 
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with complex coefficients and a multipath canceller comprising an array of complex delay-Doppler 
shift coefficients, delayers of one sample and shifters of one frequency increment, said array 
producing a residual from which the interferences from said fixed reflectors, said moving reflectors 
and said targets have been removed, said multipath canceller comprising a minimizer minimizing 
the mean square of said residual over the said whitening filter coefficients and the said multipath 
canceller delay-Doppler coefficients, said surveillance subsystem comprising, for each receiver, a 
target measurer giving differential range and differential range rate from said delay-Doppler 
coefficients, said target measurer, if more than one antenna is provided, giving target angle from the 
relative phases of said complex delay-Doppler coefficients, said surveillance subsystem comprising 
a target tracker utilizing the said measured differential range, said differential range rate, said angle 
and the position and velocity of the platform on which said antennas are mounted, said position and 
said velocity being determined by a navigation subsystem, said tracker locating and tracking said 
target or targets. 

12. A surveillance subsystem for use in a surveillance system with targets which are moving 
reflectors, with interfering fixed reflectors and with interfering moving reflectors, said system 
comprising one or more stationary or moving transmitters, one or more stationary or moving 
antennas and one or more receivers, said surveillance subsystem comprising a multipath reduction 
subsystem for each receiver, said multipath reduction subsystem comprising an analog to digital 
converter, a real to complex converter, a whitening filter with complex coefficients and a multipath 
canceller comprising an array of complex delay-Doppler shift coefficients, delayers of one sample 
and shifters of one frequency increment, said array producing a residual from which the 
interferences from said fixed reflectors, said moving reflectors and said targets have been removed, 
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said multipath canceller comprising a minimizer minimizing the mean square of said residual over 
the said whitening filter coefficients and the said multipath canceller delay-Doppler coefficients, 
said surveillance subsystem comprising, for each receiver, a target measurer giving differential 
range and differential range rate from said delay-Doppler coefficients, said target measurer, if more 
than one antenna is provided, giving target angle from the relative phases of said complex delay- 
Doppler coefficients, said surveillance subsystem comprising a target tracker utilizing the said 
measured differential range, said differential range rate, said angle and the position and velocity of 
the platform on which said antennas are mounted, said position and said velocity being determined 
by a navigation subsystem, said tracker locating and tracking said target or targets. 

13. A surveillance method for use in a surveillance system with targets which are moving reflectors, 
with interfering fixed reflectors and with interfering moving reflectors, said system comprising one 
or more stationary or moving transmitters, one or more stationary or moving antennas, one or more 
receivers each with an analog to digital converter, said surveillance method comprising reducing 
multipath for each receiver, by converting real sampled data to complex sampled data, whitening 
said complex sampled data by filtering with complex coefficients and cancelling said multipath by 
delaying and frequency shifting the residual under control of the complex delay-Doppler shift 
coefficients producing said residual from which the interferences from said fixed reflectors, said 
moving reflectors and said targets have been removed, said whitening filtering coefficients and said 
multipath cancelling delay-Doppler shift coefficients being determined by minimizing the mean 
square of said residual, said surveillance method comprising, for each receiver, measuring said 
target or targets by obtaining differential range and differential range rate from said delay-Doppler 
shift coefficients, and if more than one antenna is provided, obtaining target angle from the relative 
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phases of said complex delay-Doppler shift coefficients, said surveillance method comprising 
locating and tracking said target or targets utilizing the said measured differential range, said 
differential range rate, said angle and the position and velocity of the platform on which said 
antennas are mounted, said position and said velocity being determined by a navigation subsystem. 
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